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Abstract 

We propose a form for the action of a relativistic particle subject to a positional force that is invariant under 
time reparametrization and therefore allows for a consistent Hamiltonian formulation of the dynamics. This 
approach can be useful in the study of phenomenological models. Also the Dirac and Klein-Gordon equation 
differ from the standard formulation, with corrections of order {E — m)/m in the energy spectra. 
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1. Introduction 

Analytical mechanics is at the basis of modern theoretical physics. In particular, Hamiltonian methods 
are largely employed in several applications and are essential for the formulation of the quantum theory. 
However, not all systems have been consistently described in Hamiltonian terms. For example, if one 
excludes fundamental interaction, like Maxwell theory, whose coupling to particles is dictated by gauge 
invariance, in special relativity it is not easy to introduce particle interactions that preserve the Lorentz 
and the reparametrization invariance of the action. While for phenomenological interactions the request 
of Lorentz invariance is not compelling, the breakdown of the reparametrization invariance has serious 
consequences on the possibility of defining a consistent Hamiltonian for the model under study. The lack 
of a Hamiltonian formulation also implies that in quantum theory the Klein-Gordon or Dirac equations for 
generic interactions have to be postulated rather than derived from the correspondence principle. 

These problems arise especially for models like that of an external central force acting on a particle, when 
the force is fixed a priori and is not determined by field equations derived from a variational principle. An 
important case is that of the harmonic oscillator, whose relativistic formulation is problematic [1]. Usually, 
in such cases the Hamiltonian is defined only for a specific choice of the time coordinate. 

In this paper, we propose a different solution that, slightly modifying the coupling, makes the action 
reparametrization invariant. This allows to obtain a Hamiltonian by means of the usual Dirac formalism 
for constrained systems. The corrections with respect to the standard formalism are of order {E — m)/m, 
where E is the relativistic energy and m the mass of the particle. Moreover, a really covariant Klein-Gordon 
or Dirac equation can be obtained. In the case of the harmonic oscillator, a proposal leading to equivalent 
results was advanced long ago in [1], starting from different considerations. 

The paper is organized as follows: in sect. 2 we review the Lagrangian and Hamiltonian formalism for 
a free particle. In sect. 3 we describe our proposal and compare it with the standard formulation of the 
relativistic interacting particle. In sect. 4 we give some elementary instances of application of the formalism. 
In sect. 5 the changes in the Klein-Gordon equations are discussed and illustrated with simple examples in 
sect. 6. 



2. Free particle 

In special relativity, the motion of a free particle in coordinates is determined by the variation of the 
action 

S = j CdT = -mj ^/{±pfdT = -mj ds, (1) 

where C is the Lagrangian density, t an arbitrary evolution parameter, a dot denotes a derivative with 
respect to r, m is the rest mass and s the proper time on the trajectory. The action is invariant under 
Lorentz transformations and under reparametrizations r t'{t). The Euler-Lagrange equations for the 
action (1) read 

d 

m ; = 0. (2) 

dT 

The dynamics can also be written in Hamiltonian form [2]. The momentum conjugated to is defined 

as 

dC da;'* 
Pu = TTT- = m — , = m——, (3) 
dxf' V^i^ ds ^ ' 

and satisfies the constraint 

pI = m\ (4) 

Since the action is reparametrization invariant, the Hamiltonian % vanishes, and the action can be 
written in first-order form as 



s = j drix^'pf, -xn), n = p' 



2 ™2 



m\ (5) 



where A is a Lagrange multiplier enforcing the constraint (4). 
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The Hamilton equations ensuing from the action (5) are 



i/x = {Xfj., XH} = 2Xp^, = {p^, XH} = 0, (6) 

from which, comparing with (3), follows that 

A=i^^ (7) 

In Dirac formalism, the constraint (4) is first order. One can therefore reduce the system by choosing 
a gauge, i.e. fixing the time coordinate. 

The standard choice is t = x°, namely one identifies the evolution time t with the coordinate time xq. 
Since {t, XH} = 2Xpo, no secondary constraints arise and one easily checks that 



/ Pi Pi 
X, = 



^:=o, a=-L0t;^ 1 (8) 



Po ^pi + m? ' 2m V « 2po 

where i = 1, 2, 3 and a prime denotes a derivative with respect to t. The reduced action is then 



S = J dt ^x'^Pi 



Pl + m?\, (9) 



and the 3-dimensional effective Hamiltonian H = po = + m? can be identified with the energy of the 
particle in the laboratory frame. 

Other choices of the gauge are however possible. For example, with the choice t = x^p^, the evolution 
time coincides with the proper time of the particle. 

3. Particle in an external potential 

The addition of an external potential acting on the free particle presents some problems and to our 
knowledge has not been discussed in depth. In the standard formalism, one simply adds to the action a 
potential term 

Sint = - j V{x,x)dT. (10) 

In general the potential breaks the Lorentz invariance. Moreover, unless V is homogeneous of degree one in 
the velocity, as in the Maxwell case when V = i;^A^(a;), Sint is not invariant under reparametrization, and 
the Dirac formalism cannot be consistently applied for writing down the equations of motion in Hamiltonian 
form. 

In the following we discuss the case of positional potentials V = V{x) that do not depend on the velocity 
of the particle. These usually arise as phenomenological potentials, like those associated to elastic forces. 
For positional potentials, the interaction term (10) leads to the standard equations of motion 

A. jfiL - d^v = 0. (11) 

Multiplying eq. (11) by x^ and integrating, one can show that the total energy 

E = + V (12) 

is conserved. Since V = V{x), the momenta are given by (3) as for the free particle and satisfy the constraint 
(4). If one tries to apply the Dirac formalism to the action so defined, one obtains that V must be a constant. 
This is a consequence of the lack of reparametrization invariance of (10). 
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In order to avoid these problems and to preserve reparametrization invariance, we propose an interaction 
term of the form 



Sint = - I ^Ai^y{^)dT, (13) 

leading to the total action 



S=- j sJixPf (m + V) dr. 
The Euler-Lagrange equations for the action (14) read then 

dr ^/p^ ^ 
If V does not depend on i;'', cq. (15) can also be written 



(14) 



(15) 



(m + y)A_£^ + -^---^-^^-^-(f)^^-^=0. (16) 
^ 'dr v^(i^ ^/(i^ 

One can again find the energy integral by multiplying eq. (15) by and integrating, obtaining 

E=^. (17) 
To pass to the Hamiltonian formulation, we compute the momentum conjugate to x^, 

which is subject to the constraint 

Pl = {m + Vf. (19) 
It follows that the action can be written in Hamiltonian form as 

5 = j drix^'pf,- Xn), n=pl - {m + Vf (20) 

where A is a Lagrange multiplier enforcing the constraint (19). 
The Hamilton equations read 

x^ = {x^, XH} = 2Xp^, = {p^, Xn} = -2A(m + V)d^V, (21) 

from which follows that 

^~2im + Vy ^^^^ 

In our formalism, the potential essentially plays the role of a variable mass added to the rest mass of the 
particle. This is in accord with special relativity, since the external potential is fixed and does not change 
with the motion of the particle. So its field cannot contribute to the energy balance in any other way than 
effectively modifying the mass of the particle. 

Again, one can choose a gauge in order to reduce to a three-dimensional problem with external time. 
Taking t = x°, one has {t — xo, Xy,} = 2XpQ, and hence no secondary constraints arise if the potential does 
not depend on xo- Moreover 

p', = -ni±^d.V, A= y^ =-L. (23) 



' Po ^pf + {m + V)^' ' Po ' ' 2{m + V) 2po' 
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In this case the reduced action is 



dt 



(x',pi - ^pj + {m + Vy^ . (24) 



The 3-dimensional Hamiltonian H = p^ + (m + VY" represents the energy in the laboratory frame. It 
should be compared with the expression adopted in the standard formalism, namely H = ^/pf+rn^ + V. 
This is obtained from the reduced lagrangian of the standard formalism in the gauge t = xq and has no 
covariant meaning. The two Hamiltonians differ by terms of order {E — m)/ m. 

For m 3> p,;, rn >■ V, the expansion of both Hamiltonians gives rise to the classical expression plus 
higher-order corrections, 

E^m+-^ + V + ... (25) 
Zm 

However, with our Hamiltonian higher-order corrections are present in the expansion also for the potential 
term. 

Like for the free particle, also in presence of a potential term, one can choose a different gauge. For 
example, in 1+1 dimensions, the evolution time can be made to coincide with the proper time of the particle 

for the choice t = {m + V)(l + f P- 
4. Examples 

We give two elementary examples that we solve in the Lagrangian formalism, and compare them with 
the results obtained using the standard formalism. The potentials are not Lorentz invariant and therefore 
the results hold only in the rest reference frame. However the equations of motion are reparametrization 
invariant. 

4-1 Harmonic oscillator 

For a 1-dimensional system, the equations of motion are immediately integrated. Choosing t = xq, the 
conserved energy (17) can be written 

^=^^- (26) 



VT- X' 



Easy calculations lead to the differential equation 

1 



x' = - ^E^ - (m + y)2. (27) 
E 

that can be immediately integrated. In the case of a harmonic oscillator, V = ^kx^, the same equation has 
been obtained in [1], starting from a different perspective. 
The solution of (27) is given by [1] 




with period 



^^sdU/^t,::^), (28) 



^ , E/E-m\ „ /m/ b E - m , 



where sd and K are elliptic functions. Contrary to nonrelativistic mechanics, the period of the oscillations is 
not constant, but depends on the energy. We have written down the first term of its expansion in {E — m)/m 
around the nonrelativistic value. 

Note that in the standard formalism, the equation of motion would have read 



,30) 
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Also this equation can be solved in terms of elliptic functions [3] . In particular, the period of the solution is 



whose dependence on energy differs from (29). 

4-2 Kepler problem 

It is well known that the Kepler problem can be reduced to a one-dimensional problem because of the 
presence of two integrals of the motion related to the conservation of the angular momentum. This implies 
that the motion occurs on a plane and that the component of the angular momentum orthogonal to the 
plane is conserved. Taking t = xq and polar coordinates r and ip on the plane of the motion, from (17) one 
has 

E= , "^ + ^ (32) 
while from the conservation of the norm of the angular momentum follows that 



is a constant. After standard computations from (32) and (33) one obtains the equation of the orbits. 
Defining u=- and V = —— = —au, eq. (34) becomes 



2 

= ^ [E^ - (to - auf] (35) 



dip J 



whose solution is given by 



with 



r = - = -. T' (36) 

u 1 + e cos q((p — (fo) 



, Q!2 P + a'^ E r (£2 _ ^2) 12 



P ' am ' TO V a^E^ 

These values should be compared with those obtained using the standard formalism, namely 



a2 P - q;2 m l (E'^ — m'^)P 

In the formulae above, e is the eccentricity of the orbit, while the parameter p is related to the size of the 
orbit and q to the perihelion advance S = 27r(g~^ — 1). The angle S has opposite value in our formalism 
with respect to the standard one. Of course, the correct value is that given by general relativity and is three 
times the standard one. The other parameters of the orbit calculated in our formalism differ as usual for 
terms of order {E — to) /to from the standard ones. 

5. Quantization 

A first quantized relativistic equation can be obtained as usual starting from the Hamiltonian (20). The 
Klein-Gordon equation is obtained with the substitution ihdij,, namely 

<S^-S)h-(---)V = 0. (3.) 
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Setting ij; = (j)e '-^t/^, it reduces to 

-n^^ + [{m + Vf-E^]<t> = Q, (40) 

that has the form of a Schrodinger equation with -Be// = E"^ — m^, Vgff = 2mV + V'^. Note that if the 
potential V is positive definite, also Vgff is. This is not true in the standard formalism, leading to problems 
with stability. 

In fact, in the standard formalism, the Klein-Gordon equation is assumed to be 

-n'^ + [m'-{E-Vr]cP = 0. (41) 

It also has the form of a Schrodinger equation with -Be// = E"^ — m^, Vgff = 2EV — V"^. Note that this 
equation cannot be derived from a classical Hamiltonian, since that is not well defined. Moreover, the 
presence of E in the effective potential makes it difficult to find the spectrum of the energy explicitly. 
Finally, we notice that also the Dirac equation can be modified in accord with (39), as 

[ih^^d^, - (m + y)]V = 0. (42) 



6. Examples 

We give two simple examples of solution of the Klein-Gordon equation. In the first case, only a per- 
turbative solution is possible, while the second is constructed so that an exact solution can be found. For 
simplicity, in this section we set ?i = 1. 

6.1. Quantum harmonic oscillator 

The Klein-Gordon equation (39) for a one-dimensional harmonic oscillator with potential V = ^kx^ 
reads ^ 

^ = {jjLx'^ +mkx'^ - e)(j), (43) 
ax-'' 

where e = E^ — m'^, /z = fc^/4, and has the form of a nonrelativistic biquadratic oscillator. Its energy 

spectrum can be obtained using standard perturbation theory, taking the quartic term as a perturbation of 
the nonrelativistic harmonic oscillator. The calculation has been performed in [1] and gives at first order 



E„=± 



m + 



-Min +n+- 



(44) 



where ri is an integer. 

In the standard formalism, the Klein-Gordon equation (38) reads instead 



dx'^ 



= {-^x^ + Ekx^ - e) 



(45) 



with e and /i as before, and the effective potential is no longer positive definite, leading to possible instabilities 
[4]. Neglecting this problem, one can obtain as before the energy spectrum, which at first order in Ti, and 
neglecting corrections of order {E — m)/m, yields the opposite sign for the correction, namely 



En 



± 



/in +n + 



1 



(46) 



7 



6.2. Exactly solvable potential 

An instance in which a solution of the Klein-Gordon equation (39) can be found explicitly is 



V = m 



1 



cos a: 



1 . 



(47) 



This is a potential well whose value is zero at the origin and infinite at a; = The Klein-Gordon equation 
reads 



m 



Defining z = sin x, eq. (48) can be put in the form 



dz'^ dz 



= 0. 



E' 



= 0, 



.(1-^2)2 i_^2_ 

that in turn can be reduced to a standard hypergeometric equation with solution 

4> = cos'^x F{u + E,i/- E,l'+-,- 



m 



(49) 



(50) 



where F is a hypergeometric function and v = 5 (l + VY+'ImF). Imposing the vanishing oi (j) at x = 
one finds the eigenfunctions 

0„ = cos'^a; C^'^^(sina;), (51) 



with Cn^ Chebyshev polynomials, and the energy spectrum 



E = ± 



m 



1/2 



(52) 



where n is an integer. 

In the standard formalism, the Klein-Gordon equation cannot be solved analytically and therefore we 
do not pursue its investigation here. 



7. Conclusions 

We have shown that it is possible to formulate the problem of the motion of a particle in a positional 
potential in special relativity in such a way to preserve reparametrization invariance, and hence to give a 
consistent Hamiltonian formulation in terms of the Dirac formalism for constrained systems. Of course, 
potentials of the kind investigated in this paper do not occur in fundamental interactions, whose action is 
already invariant under reparametrization, but can be of interest for phenomenological models. 
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